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ARSTRACT 

Sonic  full  waveforms  were  recorded  in  a  shale- 
limestone  interval  with  an  experimental  12-channel 
acoustic  logging  device.  Well-defined  shear  (pseudo- 
Rayleigh)  waves  were  observed  throughout  most  of  the 
142  m  interval  and  were  used  to  distinguish  lithologic 
boundaries  and  zones  ofVracturing  as  interpreted  by 
acoustic  borehole  televiewer  and  standard  geophysical 
logs.  The  high  signal-to-noise  ratio  of  the  waveforms 
permitted  meaningful  shear  velocities  to  be  obtained  by 
12-fold  semblance  calculations  in  all  but  the  washed-out 
intervals  in  this  well.  A  correlation  of  changes  in  the 
V-  l's  ratio  with  depth  aided  in  identifying  some  compo¬ 
sitional  changes  in  the  well,  but  such  a  correlation  could 
not  be  used  to  distinguish  fracturing  in  the  limestone. 
Shear-wave  amplitudes,  however,  were  effective  for  iden¬ 
tifying  lithologic  changes,  as  well  as  fractures,  from  the 
sonic  log.  Shear  amplitudes  were  lower  in  the  shale  than 
in  the  limestone,  and  shear-wave  attenuation  increased 
in  fractured  zones.  There  was  no  strong  correlation  be¬ 
tween  the  degree  of  fracturing  and  the  attenuation  of 
Stoneley  waves.  In  general,  shear-wave  processing  of 
sonic  logs  is  recommended  for  interpretation  of  lithol¬ 
ogy  and  fracturing  in  any  well  in  which  full  waveforms 
arc  recorder  and  shear  waves  propagate. 


INTRODUCTION 

In  theoretical  studies,  the  V  j Vt  ratio  has  been  related  to 
lithology,  crack  density,  and  pore-fluid  saturation  (O'Connell 
and  Budiansky,  1974;  Toksoz  et  al.„  1976).  These  models  gen¬ 
erally  predict  that  V  /  Vs  will  increase  as  the  spatial  crack  den¬ 
sity  increases  in  a  homogeneous  medium,  where  the  cracks  are 
much  smaller  than  the  seismic  wavelength.  Laboratory  investi¬ 
gations  have  also  convincingly  shown  that  the  presence  of  flat 
cracks  or  pores  in  a  homogeneous  medium  significantly  re¬ 


duces  its  elastic  constants,  and  hence  its  seimsic  velocities 
(Birch.  I960;  Nur  and  Simmons.  1969;  Brace  et  al..  1972; 
kuster  and  Toksoz.  1974;  Mavko  and  Nur.  1979).  Tatham 
(1982)  suggests  that  the  dependence  of  V  V\  on  crack  density 
is  even  stronger  than  on  the  elastic  constants  of  the  minerals 
comprising  the  matrix,  but  an  association  of  crack  density 
with  lithology  may  enable  a  correlation  between  V  iVs  and 
lithology.  On  the  other  hand,  Wilkens  et  al.  (1984)  have 
shown,  in  laboratory  measurements  on  limestone  samples,  a 
stronger  dependence  of  V  1',  on  carbonate  content  than  on 
crack  closure.  Although  the  relative  contributions  of  compo¬ 
sition  and  microcracking  on  elastic  properties  are  still  uncer¬ 
tain.  empirical  relationships  have  commonly  been  applied  to 
sonic  logging  data  to  interpret  lithologic  changes  from  Vp  f' 
(Pickett.  1 963 ;  Nations.  1974;  Kithias.  1976;  Eastwood  and 
Castagna.  1982).  These  applications  are  limited,  however,  be¬ 
cause  no  unique  relationships  have  been  found  (Moos  and 
Zoback.  1983). 

In  this  paper,  (tie  sonic  waveform-derived  Vp  Ij  ratio  is  cor¬ 
related  with  lithology  and  fracturing  in  a  shale-limestone  se¬ 
quence.  Considering  the  likelihood  of  nonunique  results,  we 
also  investigate  the  amplitude  and  frequency  of  borehole  wave 
modes  as  simple  and  sensitive  indicators  of  changes  in  the 
borehole  environment.  Paillet  and  White  (1982)  fully  discuss 
the  theory  of  an  acoustic  wave's  sensitivity  to  the  borehole 
environment.  In  particular,  relative  changes  in  shear-wave  am¬ 
plitudes  are  largely  dependent  upon  guided  (pseudo-Rayleigh) 
waves  in  the  borehole,  and.  as  such,  are  affected  by  environ¬ 
mental  features  which  are  difficult  to  identify  from  tradition¬ 
ally  processed  sonic  logs.  The  processing  described  in  this 
paper  is  effective  for  the  full-waveform  data  set  of  any  sonic 
log  in  which  shear  waves  propagate. 

DATA 

In  this  study.  Schlumbcrger's  experimental  12-channel  sonic 
sonde  was  used  to  record  acoustic  waveforms  from  source 
depths  of  1469.4  to  1327.4  m  in  a  limestone-shale  sequence  in 
the  Anadarko  Basin,  Oklahoma.  The  source  is  near  the 
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bottom  of  the  sonde  and  3.96  in  below  the  nearest  of  a  linear  Core  descriptions  in  a  nearby  well  support  the  identification 

array  of  12  equally  spaced  receivers  spanning  1.68  m  The  of  the  shale,  limestone,  and  calcareous  shale  intervals  from  the 

sonde  was  pulled  upward  at  a  constant  rate  from  the  bottom  logs.  The  shales  are  highly  interbedded  and  brecciated.  Chert 

of  the  hole,  and  a  set  of  waveforms  was  recorded  for  each  nodules  are  observed  in  layers  up  to  15  cm  thick  in  the  upper 

source  location  at  1 5  cm  depth  intervals.  The  source  generates  10  m  of  the  limestone,  and  extensive  fracturing  is  observed 

an  acoustic  signal  having  a  peak  energy  level  in  the  5  to  15  throughout.  These  features  are  evident  in  an  acoustic  tele- 

kHz  frequency  band  (Kimball  and  Marzetta,  1984).  The  viewer  log  in  the  limestone  interval,  where  fractures  intersect- 

source  was  fired  once  for  each  receiver  recording.  Rapid  firing  ing  the  borehole  appear  as  low-reflectance  (dark)  bands  in  the 

of  the  source  for  the  succession  of  12  receivers  resulted  in  a  image  (Figure  2).  Near  1433.0  m,  the  observed  fractures  are 
negligible  change  in  source  location  during  recording.  Each  presumably  the  largest  in  the  logged  interval  and  have  a  pre¬ 
waveform  was  digitized  uphole  at  a  10  ps  sampling  rate.  dominantly  vertical  orientation.  The  televiewer  traveltime  log. 

The  compressional,  shear,  and  Stoneley  wave  mode  arrivals  also  shown  in  Figure  2,  corresponds  to  the  relative  position  of 

are  prominent  and  are  identified  for  waveforms  recorded  at  the  televiewer  to  the  borehole  wall.  The  traveltime  log  shows 

the  near  receiver  in  Figure  1.  The  near-receiver  waveforms  the  tool  is  slightly  off  center,  but  there  is  little  ellipticity  of  the 

have  been  compensated  for  low-frequency  instrument  biases  cylindrical  hole,  so  that  sonic  waveform  distortion  caused  by 

and  gain  controlled  to  equalize  the  amplitude  of  late-arriving  the  ellipticity  of  the  borehole  is  negligible, 

modes.  Other  geophysical  logs  recorded  in  this  well  are  also 

displayed  in  Figure  1.  The  interpretation  of  these  logs  delin-  DATA  PROCESSING 

eates  a  shale  interval  above  1393.3  m,  a  limestone  interval 
below  the  shale  to  1448.2  m,  and  a  calcareous  shale  to  the  Velocity  analysis 

bottom  of  the  well.  Shale  intervals  generally  have  low  induc¬ 
tion  and  high  gamma-ray  readings,  and  in  certain  zones  high  In  the  past,  semblance  has  been  applied  extensively  to  seis- 

caliper  readings  indicate  the  presence  of  washouts.  In  the  lime-  mic  reflection  data  to  obtain  stacking  velocities.  By  adapting 

stone,  neutron  and  density  porosity  measurements  are  rela-  to  the  linear  moveout  of  refraction  data,  multichannel  some 

tively  consistent,  indicating  good  hole  conditions,  and  the  cali-  waveforms  can  be  analyzed  using  semblance  to  obtain  interval 

per  log  in  this  interval  is  constant,  but  20  percent  lower  than  velocities.  The  semblance  Su  was  described  by  Taner  and 

in  the  shale.  Koehler  (1969)  as  a  measure  of  the  coherence  of  multichannel 
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Fig.  1.  Geophysical  logs  from  1327.4  to  1469.4  m  depth.  Gain-adjusted  near-receiver  waveforms,  compressional  and 
shear  velocity  logs,  and  V JVt\  induction,  gamma-ray,  and  caliper  logs;  and  density  (solid)  and  neutron  (dashed) 
porosity  curves  Three  washouts  of  the  borehole  occur  in  the  shale  above  139t  t  m  From  core  analysis,  chert  iayers 
weie  identified  in  the  upper  10  m  of  the  limestone  and  extensive  fracturing  occurs  throughout.  A  porous  calcareous 
shale  underlies  the  limestone. 
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REFLECTANCE  TIME  REFLECTANCE  TIME 


Fig.  2.  Borehole  televiewer  reflectance  and  traveltime 
images  in  the  limestone.  Dark  vertical  banding  indicates  low 
reflectance  due  to  intersecting  fractures.  The  most  intense  frac¬ 
turing  occurs  near  1433  0  m. 
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Fig.  3.  Amplitude  spectra  of  isolated  compressional.  shear, 
and  Stonelev  arrived  now  a  w*<eform  recorded  at  1469.4  m. 
Stoneley  energy  has  significantly  lower  peak  frequency  than 
does  either  the  compressional  or  shear  energy. 


signals,  normalized  by  their  autocorrelations  at  each  time  lag, 
and  given  by 


M  =  the  number  of  traces, 
n  =  1/2  N, 

N  =  the  number  of  samples, 
k  —  the  center  sample  in  the  time  window. 


and 


J]j  =  the  ith  sample  of  the  /th  trace. 

Semblance  ranges  between  0  and  1  for  random  and  perfectly 
coherent  signals,  respectively,  and  can  be  used  to  calculate 
velocity  by  identifying  the  time  lag  at  which  the  peak  sem¬ 
blance  value  occurs.  The  time  lag  equals  the  receiver  offset 
divided  by  the  propagation  velocity;  the  peak  semblance  value 
is  directly  related  to  the  reliability  of  the  calculation. 

To  enhance  the  accuracy  of  the  semblance  velocity  measure¬ 
ments,  time  windowing,  frequency  filtering,  and  velocity  filter¬ 
ing  were  used  to  remove  interfering  signals  and  noise  for  each 
individual  wave  mode.  Figure  3  shows  the  bandwidths  of  am¬ 
plitude  spectra  of  time-windowed  compressional.  shear,  and 
Stoneley  arrivals.  Since  shear  and  compressional  bands  over¬ 
lap.  frequency  filtering  is  not  useful,  but  time  windowing  is 
generally  adequate  to  isolate  arrivals.  The  Stoneley  spectral 
band  is  sufficiently  different  that  band-pass  filtering  alone  can 
separate  Stoneley  waves  from  interfering  arrivals.  The  results 
of  filter  combinations  are  shown  in  Figure  4.  A  suite  of  gain- 
compensated  waveforms  recorded  at  a  source  depth  of  1405.0 
m  are  displayed  on  the  left-hand  side  of  the  figure,  and  com¬ 
pressional.  shear,  and  Stoneley  arrivals,  as  identified  by 
arrows,  have  progressively  increasing  moveout.  The  second 
suite  from  the  left  shows  the  application  of  a  band-pass  filter 
to  discriminate  against  high-frequency  and  low-frequency 
noise  and  Stoneley-wave  energy.  A  velocity  filter  is  applied  in 
the  third  suite  to  discriminate  against  high-velocity  compres¬ 
sional  energy.  In  the  final  suite,  the  combination  of  the  band¬ 
pass  and  velocity  filters  used  previously  isolates  low -velocity 
shear  energy.  Other  combinations  of  filters  can  be  designed  to 
remove  different  signals  or  to  isolate  weak  arrivals. 

After  processing  to  enhance  the  waveforms,  compressional, 
shear,  and  Stoneley  arrivals  are  isolated  by  picking  the  ap¬ 
proximate  start  time  of  each  arrival  and  using  window  lengths 
of  approximately  1.5  periods.  Semblance  is  then  calculated 
and  normalized  to  its  maximum  value  High-quality  velocities 
with  relatively  high  semblance  were  obtained  over  most  of  the 
logged  interval  for  all  three  wave  modes.  The  velocities  and 
normalized  semblance  values  are  displayed  for  each  mode  in 
Figure  5.  Summing  of  the  semblance  curves  between  adjacent 
depth  points  improves  the  stability  of  the  velocity  calculation, 
with  some  loss  of  resolution.  However,  the  velocity  resolution 
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Fig.  4.  Processed  waveforms  across  the  receiver  array  at  a  source  depth  of  1405.0  m.  Waveforms  are  first  gain- 
controlled  and  bias-compensated,  then  band-pass  filtered  to  remove  high-frequency  and  low-frequency  noise  and  to 
separate  the  Stoneley  energy,  velocity-filtered  to  remove  compressional  energy,  and  finally  band-pass  and  velocity 
filtered  to  isolate  shear  energy. 


for  thin  beds  is  not  significantly  affected  if  the  summing  oper¬ 
ator  is  shorter  than  the  receiver  span.  For  these  displays,  the 
semblance  is  summed  over  0.76  m  intervals.  The  resulting 
semblance  is  low  for  all  wave  modes  only  in  the  washed-out 
intervals  between  1390.2  and  1375.0  m,  1359.8  and  1347.6  m, 
and  1333.8  and  1327.4  m  (Figure  5).  In  those  zones,  fluid, 
shear,  and  Stoneley  velocities  are  all  near  1500  m/s,  so  that 
their  arrivals  interfere  with  each  other  and  result  in  low  sem¬ 
blance.  Frequency  and  velocity  filtering  to  separate  these  ar¬ 
rivals  did  not  effectively  increase  the  quality  of  the  calculated 
velocities.  The  irregular  shape  of  the  borehole  in  these  zones 
may  be  a  major  contributing  factor  to  the  difficulty  in  separ¬ 
ating  wave  modes. 

Figure  6  illustrates  waveforms  recorded  across  a  sharp  shale- 
limestone  boundary.  Low-amplitude  compressional  waves  are 
not  visible  to  the  eye  at  this  plotting  scale.  Shear-wav?  move- 
out  dramatically  decreases  as  the  receiver  array  progresses 
across  the  velocity  contrast  and  from  left  to  right  across  the 
figure.  Shear  velocity  increases  by  about  1000  m/s  across  the 
boundary,  and  because  the  velocity  difference  is  large,  a  sem¬ 
blance  peak  is  calculated  corresponding  to  tb-  velocity  on 
each  side.  The  layer  having  more  receivers  will  have  higher 
peak  semblance.  In  the  case  of  narrower  velocity  differences 
between  layers,  semblance  correlation  results  in  averaging  of 
the  velocities  over  the  receiver  span. 

Amplitude  analysis 

Amplitudes  of  the  waveforms  were  calculated  by  two  differ¬ 
ent  methods.  In  the  first  method,  a  time-domain  measure  of 


amplitude,  the  amplitude  envelope  was  calculated  using 
=  f/2(t)  +  /«(')] 1  2-  where  /„(()  equals  the  Hilbert  trans¬ 
form  of  /(t)  (Bracewell.  1965).  The  computed  amplitude  en¬ 
velope  is  shown  in  color  in  Figure  7  for  near-receiver  and 
far-receiver  waveforms.  The  vertical  axis  corresponds  to  depth 
in  the  well,  and  the  1 .68  m  offset  between  near  and  far  receiv¬ 
ers  is  not  aligned  in  order  to  maintain  a  common  source 
location  for  each  depth  point.  The  time  scale  is  truncated  to 
exclude  the  high-amplitude  Stoneley  arrival  from  scaling. 
Compressional  and  shear  attenuation  is  apparent  by  a  loss  of 
amplitude  between  near  and  far  receivers.  Attenuation  is 
greatest  for  shear  waves  near  1433.0  m.  where  losses  up  to  70 
percent  occur  across  the  array,  although  elsewhere,  shear  at¬ 
tenuation  is  about  50  percent  across  the  array.  The  variation 
of  the  amplitude  envelope  over  the  interval  suggests  a  corre¬ 
lation  between  attenuation  and  lithology,  fracturing,  and  bore¬ 
hole  conditions. 

The  second  measure  of  amplitude  is  derived  from  the  cross¬ 
power  spectrum  of  the  near-receiver  and  far-receiver  wave¬ 
forms.  The  crosspower  spectrum  can  be  calculated  by  the 
Fourier  transform  of  the  ci osscorrelation  of  time  signals  (Rob¬ 
inson  and  Treitel,  1980).  The  crosscorrelation  of  signals  x(t) 
and  y(t)  is  defined  by 

I  T 

rxy(k)  =  lim  —  -  £  *<'  +  *>><»)- 

r—  x  ^ ‘  “*■  *  *=  - r 

where  T  is  the  number  of  time  samples,  and  k  is  the  time  lag. 
The  Fourier  transform  of  r^fk)  may  be  expressed  as 

R„(w)  =  |  rxy(k]e~l°*  dk. 
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Flo.  5.  Calculated  velocities,  normalized  semblance,  and  time 
picks  of  isolated  compressional,  shear,  and  Stoneley  wave¬ 
mode  arrivals.  Semblance  is  summed  over  0.76  m  intervals  to 
improve  the  stability  of  the  velocity  calculation.  Low  sem¬ 
blance  occurs  in  washouts  and  at  severe  velocity  contrasts. 


Peak  magnitude  and  frequency  are  calculated  from  Rxy  to)  in 
order  to  correlate  better  sonic  amplitude  and  frequency  data 
with  lithology,  fracturing,  and  borehole  conditions.  These  re¬ 
sults  are  discussed  below. 

RESULTS 

The  yr/y,  ratio 

In  Figure  8,  the  results  of  compressional  and  shear  velocity 
analyses  have  been  crossplotted  and  grouped  for  detailed 
lithologic  interpretation  by  a  depth  correlation  between  the 
crossplot  and  the  logs.  The  slope  of  the  data  in  the  crossplot 
represents  the  Vp/Vs  ratio  and  varies  between  about  1.50  and 
2.0  over  the  interval.  Subgroupings  of  the  data  labeled  A 
through  J  illustrate,  for  example,  that  washout  intervals  (J) 
have  constant  shear  velocities  equal  to  the  fluid  velocity  (1500 
m/s),  and  therefore  data  from  those  intervals  are  believed  to  be 
unreliable.  Compressional  and  shear  velocities  are  greater  by 
about  1000  m/s  in  the  limestone  (F,  G,  H,  and  I)  than  in  the 
shale  (B.  C,  and  E),  although  the  Vp/VM  value  is  lower  in  the 
shale.  The  fractured  (G)  and  chert  interbedding  (F)  intervals 
have  lower  velocities  than  the  surrounding  limestone,  al¬ 
though  I '  IV,  remains  relatively  unchanged.  In  addition,  the 
V JVa  ratio  is  not  constant  within  the  shale  (B  and  C)  and 
calcareous  shale  (A  and  D)  intervals.  Hence,  the  Vp/V,  ratio  is 
effective  for  distinguishing  major  lithologic  changes,  but  is  un¬ 
reliable  in  washouts,  and  cannot  be  used  to  discriminate  be¬ 
tween  fractured  and  unfractured  intervals  in  this  well. 

Crosspower  spectra 

The  peak  magnitudes  of  the  crosspower  spectra  from  iso¬ 
lated  compressional,  shear,  and  Stoneley  wave  modes  are  plot- 
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Fig.  6.  Unprocessed  waveforms  recorded  at  each  receiver  for  15  cm  incremental  souice  depths  from  left  to  right  across 
a  shale-limestone  boundary.  Shear  and  Stoneley  moveouts  decrease  progressively  as  more  receivers  move  into  the 
limestone.  Low-amplitude  compressional  arrivals  are  not  visible  to  the  eye. 
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Fig.  7.  Amplitude  envelope  color  display  for  near-receiver  and  far-receiver  waveforms  at  every  fourth  source  location. 
The  color  scale  illustrates  the  relative  variation  in  amplitude  across  the  receiver  array  and  with  depth. 
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Fig.  8.  Crossplot  of  Vp  and  VS1  where  data  have  been  grouped  corresponding  to  their  depth  intervals.  Shale  (B  and  C), 
calcareous  shale  (A  and  D),  and  limestone  (F,  G,  H,  and  I)  intervals  are  distinguishable.  Washouts  (J)  show  a  constant 
shear  velocity  (1500  m/s)  which  is  unreliable.  Chert  interbedding  <F)  and  fracturing  (G)  in  the  limestone  are  dis¬ 
tinguished  by  decreases  in  velocity,  although  the  Vp!Vi  ratio  remains  unchanged. 


ted  as  a  function  of  depth  in  Figure  9.  The  caliper  log,  which 
measures  the  diameter  of  the  borehole,  is  also  displayed.  The 
compressional  amplitude  is  high  in  the  shale  and  low  in  the 
limestone  and  calcareous  shale.  Note  that  large  decreases  in 
shear  amplitude  occur  in  the  fractured  interval  near  1433.0  m 
and  at  the  sharp  shale-limestone  boundary  at  1457.3  m,  where 
there  are  no  significant  changes  in  the  borehole's  diameter. 
The  decay  of  the  shear  amplitude  in  these  zones  is  clearly 
associated  with  these  features.  In  the  overlying  shale,  however, 
low  shear  amplitudes  may  be  due  to  both  poorer  shear-wave 
propagation  and  a  larger  borehole  diameter  than  in  the  lime¬ 
stone.  The  lowest  shear  amplitudes  recorded  in  the  well  occur 
in  the  shale. 

Both  shear  and  Stoneley  amplitudes  typically  decrease  in 
fractured  intervals  (Paillet,  1980;  1981).  although  in  this  well, 
only  shear  amplitudes  appear  to  be  affected.  Stoneley  ampli¬ 
tudes  decrease  in  washouts  and  at  the  sharp  shale-limestone 
boundary  and  are  relatively  unaffected  by  fracturing  or  the 
decrease  in  the  diameter  of  the  borehole  in  the  limestone.  This 
insensitivity  to  fracturing  and  borehole  size  may  be  due  to 
amplifier  saturation  during  recording  of  the  high-amplitude 
Stoneley  wave. 

The  peak  frequencies  of  isolated  compressional,  shear,  and 
Stoneley  wave  modes  are  also  plotted  in  Figure  9.  Compres¬ 
sional  and  shear  peak  frequencies  are  nearly  equal  and  de¬ 
crease  from  about  9  kHz  in  the  limestone  to  about  7  kHz  in 
the  shale.  Because  the  frequency  contents  of  compressional 
and  shear  waves  are  governed  largely  by  the  resonant  fre¬ 
quency  of  guided  waves  in  the  borehole,  any  variation  in  the 
peak  frequency  can  be  attributed  primarily  to  changes  in  the 
size  of  the  borehole  (Paillet  and  White,  1982).  Stoneley-wave 
frequencies  remain  constant  over  the  well,  which  may  also  be 


a  result  of  recording  saturation.  We  place  lower  confidence  in 
the  frequency  measurements  in  washouts,  however,  since  com¬ 
pressional,  shear,  and  Stoneley  arrivals  are  all  severely  attenu¬ 
ated. 


DISCUSSION 

In  this  well,  the  relationship  between  the  lithology  and  Vn  'Vs 
values  calculated  from  the  sonic  logs  is  not  unique,  since  Vp i Vs 
values  alone  cannot  distinguish  all  of  the  environmental  and 
lithologic  changes.  For  example,  thin  chert  layers  and  frac¬ 
tures  in  the  limestone  affect  the  velocities  in  such  a  way  that 
VpjVs  remains  unchanged  from  that  of  the  surrounding  rock, 
but  both  Vp  and  V5  individually  decrease  in  a  similar  manner. 
A  theoretical  dependence  of  Vp/V!l  on  microcracks  and  lithol¬ 
ogy  in  homogeneous  media  will  be  invalid  if  inhomogeneities 
are  present  which  have  a  spatial  distribution  similar  to  the 
sonic  wavelength,  which  may  explain  the  fact  that  Vp  kj  is 
apparently  constant  throughout  the  limestone.  Localized 
alteration  in  crystalline  rock  may  possibly  explain  increases  in 
Vp/V,  near  fractures  (Moos  and  Zoback,  1983),  which  implies 
that  little  alteration  has  occurred  in  this  well  or  that  the  effect 
of  partially  altered  or  open  fractures  compensates  for  any 
change  in  V  fVp  due  to  alteration.  As  a  result,  the  VpjV^  ratio 
and  fracturing  in  this  well  are  difficult  to  correlate.  These 
velocity  measurements  in  fractured  limestone  suggest  that 
compositional  and  structural  effects  on  elastic  properties  are 
complicated  and  cannot  be  explained  by  existing  theories  or 
simple  empirical  models. 

However,  the  amplitudes  of  sonic  waves,  which  are  domi¬ 
nated  by  guided-wave  modes,  are  more  sensitive  than  veloci- 
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Fig.  9.  Crosspower  magnitude  and  peak  frequency  of  isolated  wave  modes  and  the  borehole  caliper  log.  Sh--ar 
amplitude  is  attenuated  in  the  shale  and  in  the  fractured  zone  of  the  limestone.  Stoneley-wave  attenuation  occurs  in 
washouts,  but  shows  little  response  to  the  decreasing  borehole  radius  or  to  fracturing  in  the  limestone. 


ties  to  fracturing  and  environmental  changes  in  this  well. 
Below  1400  m.  the  apparent  shear-wave  attenuation  can  be 
attributed  to  t!.e  fractured  zone  and  the  calcareous  shale 
boundary,  where  the  borehole's  diameter  is  constant.  In  the 
shale,  shear-wave  attenuation  is  greater  than  in  the  limestone, 
but  the  varying  borehole  diameter  precludes  a  clear  distinc¬ 
tion  between  geometric  and  lithologic  effects.  Poor  transmis¬ 
sion  of  energy  across  bedding  planes  (as  typically  occurs  in 
shale),  fracturing,  local  changes  in  pore  fluid  content,  and  late- 
arriving  reflections  from  sharp  boundaries  may  all  reduce 
shear  amplitudes,  but  such  factors  are  only  separable  from 
geometric  effects  if  the  size  of  the  borehole  is  constant.  Appar¬ 
ent  shear-wave  attenuation  in  inhomogeneous  formations  may 
therefore  result  from  several  different  mechanisms.  Nonethe- 
1  less,  using  traditional  sonic  tools  and  simple  processing  of 

f  sonic  waveforms,  shear-wave  amplitudes  may  aid  in  identify¬ 

ing  lithologic  and  structural  features  at  small  expense. 

SUMMARY 

In  this  case  study,  borehole  wave-mode  velocities  were  cal¬ 
culated  with  high  accuracy  by  a  semblance  algorithm.  Wave¬ 
mode  isolation  increased  the  semblance  and  improved  the  ve¬ 


locity  calculation.  In  washouts,  shear  energy  is  difficult  to  iso¬ 
late  due  to  the  coincidence  of  shear,  fluid,  and  Stoneley  arriv¬ 
als.  and  the  calculated  shear  velocity  is  unreliable.  The  Vp/Vs 
ratio  is  greater  in  the  shale  than  in  the  limestone  and  cal¬ 
careous  shale.  Near-vertical  fracturing  in  the  limestone  does 
not  significantly  change  the  value,  though  both  Vp  and  Fs 
decrease  in  a  similar  manner,  probably  due  to  a  combination 
of  fracturing  and  alteration  effects  on  the  sonic  velocities. 

Sonic  amplitudes  were  also  calculated  and  correlated  with 
lithologic  and  environmental  changes.  The  limestone  exhibits 
lower  compressional-wave  amplitude  and  higher  shear-wave 
amplitude  than  the  shale,  and  both  compressional  and  shear 
frequencies  are  higher  in  the  limestone  than  in  the  shale,  pri¬ 
marily  due  to  a  decrease  in  the  borehole's  diameter.  Stoneley- 
wave  amplitude  is  nearly  constant  in  both  lithologies,  even 
though  there  are  significant  changes  in  the  diameter  of  the 
borehole,  and  all  the  wave  modes  are  severely  attenuated  in 
washouts.  In  the  limestone  interval,  where  the  effect  of  the 
borehole  is  constant,  the  shear  amplitude  log  is  the  most  ef¬ 
fective  tool  for  distinguishing  lithologic  boundaries  and  zones 
of  fracturing.  In  the  shale,  where  the  borehole’s  diameter  is  not 
constant,  the  variable  geometric  effects  must  be  removed  in 
order  for  the  observed  changes  in  amplitude  and  frequency 
content  to  be  attributed  uniquely  to  lithologic  variations. 
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